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Nuclear matter is considered to be inhomogeneous at subnuclear densities that are realized in 
supernova cores and neutron star crusts, and the structures of nuclear matter change from spheres 
to cylinders, slabs, cylindrical holes and spherical holes as the density increases. In this letter, 
we discuss other possible structures, that is, gyroid and double-diamond morphologies, which are 
periodic bicontinuous structures discovered in a block copolymer. Utilizing the compressible liquid 
drop model, we evaluate their surface and Coulomb energies and show that there is a chance of 
gyroid appearance near the transition point from a cylinder to a slab. This interesting analogy 
between nuclear and polymer systems is not merely qualitative. The volume fraction at the phase 
transition is also similar for the two systems. Although the five shapes listed initially have been 
long thought to be the only major constituents of so-called nuclear pasta at subnuclear densities, 
our findings imply that this may not be the case and suggest that more detailed studies on nuclear 
pasta including the gyroid phase are needed. 

PACS numbers: 21.65.-f, 26.50.+X, 26.60.Gj, 61.20.-p 
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Nature contains numerous examples of materials with 
various morphologies at different levels of hierarchy. Nu- 
clear matter is no exception. Terrestrial atomic nuclei are 
usually almost spherical. About a quarter of a century 
ago, however, it was shown that nuclei deform at subnu- 
clear densities from spheres (SP) to cylinders (C), slabs 
(S), cylindrical holes (CH) and spherical holes (SH) as 
the density increases to ~ 10 14 g/cm 3 , the density of the 
uniform nuclear matter [l|, Owing to the similarity 
of these shapes to meat balls, spaghetti, lasagna, maca- 
roni and Swiss cheese, respectively, these structures are 
called nuclear pasta. Nuclei with the pasta structures are 
important because they are thought to actually exist in- 
side the core of supernovae and the crust of neutron stars 
and have an impact on astrophysical phenomena such as 
supernova explosions, proto-neutron star cooling, pulsar 
glitches and so forth [j, Q] • Assuming that only the five 
nuclear shapes listed above exist, many detailed studies 
have been perfomed since this discovery; for reviews, see 
Refs. 

Similar structure transitions are known to occur in 
nanostructures of block copolymers and, more interest- 
ingly, several complex structures have been discovered 
experimentally 0. Some of them have periodic bicon- 
tinuous structures such as the so-called gyroid (G) and 
double-diamond (D) morphologies, and are likely to ap- 
pear between the C and S phases. In this letter, we inves- 
tigate the possible appearance of G and D morphologies 
in nuclear pasta by employing the compressible liquid 
drop model (CLDM), which is a phenomenological model 
also used in earlier studies [H, H, @] • In this approach, nu- 
clei are approximated as charged liquid drops and the 
nuclear shape is determined by minimizing the total en- 
ergy density, in which the surface tension competes with 
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FIG. 1: Bird's-eye views of unit cubes of (a) gyroid and (b) 
double-diamond, in which bicontinuous minimum surfaces are 
shown for volume fraction u = 0.35. 



the Coulomb repulsion. It is also noted that the nuclear 
shape is a function of the volume fraction of the nucleus 
in a unit cell, u, similarly to the case of polymer shapes 
Q. We then find that there is a good chance that the G 
morphology appears near the transition point from the C 
phase to the S phase. Furthermore, the volume fraction 
of nuclei at this point is u ~ 0.35, which is very close 
to the value obtained for the polymer system, suggesting 
commom underlying physical laws between the two sys- 
tems. Because the nuclear pasta phases and the gyroid 
phase in the block copolymer have separately attracted 
many researchers' interest, our findings may offer a new 
field of interest. 

Setup. — The bicontinuous structures observed in 
block copolymers are thought to have a periodic mini- 
mum surface Q, which is a stationary surface for vari- 
ations of surface area with a fixed volume fraction. Al- 
though the D morphology was initially considered to be 
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the most likely structure, the G morphology is often 
thought to be a more probable structure. In this let- 
ter, we study the double network structures of G and D 
minimum surfaces as new types of nuclear pasta. 

It is known that these bicontinuous structures can be 
closely approximated by the following level surfaces: 
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for D where (x,y,z) are the spatial coordinates and 
a is the size of the unit cube. Nucleons are assumed to 
reside in the region that satisfies \ f(x,y,z)\ > k, which 
is called the "nucleus", k is a positive parameter that 
specifies the volume fraction, u, of the nucleus in the 
unit cube, k — > corresponds to u — ► 1 and u is a 
monotonically decreasing function of k. Note that some 
neutrons may drip out of the nucleus at high densities 
in the neutron-rich case. The shapes of the "nuclei" are 
illustrated in Fig. Q] 

Equations |T]) and @ are no longer good approxima- 
tions for very small values of u (and hence the corre- 
sponding values of k): u < 0.0354 for G and u < 0.161 
for D, since the resultant configurations are not bicontin- 
uous but pinched off. Although this poses no problem for 
the following analysis, we do not consider these config- 
urations. On the other hand, the hole structures of the 
gyroid and double-diamond, for which nucleons reside in 
the region satisfying \f(x, y, z)\ < k, are taken into ac- 
count. Again the configurations are bicontinuous only 
for u < 0.965 for the gyroid hole (GH) and u < 0.839 for 
the double-diamond hole (DH). 

In our model, the volume of the unit cell is a 3 and, as 
already stated, the nucleus occupies the fraction u. The 
number density of nucleons, n m , and the proton fraction, 
x ln , inside the nucleus as well as the number density of 
dripped neutrons outside the nucleus, n° ut , are set to be 
constant. Electrons are uniformly distributed in the cell 
and their number density is ux m n ln , obtained from the 
charge neutrality. 

Analysis. — Our analysis is an improvement of the 
previous CLDM studies [H, 0; 0- As shown below, it 
is mainly based on a geometrical argument and is insen- 
sitive to nuclear interaction models. We write the total 
energy in the unit cell, W, as 



W = W b + W s + W c , 



(3) 



where Wb, W s and Wc are the bulk energy both inside 
and outside the nucleus including the kinetic energy of 



electrons, the surface energy and the Coulomb energy, 
respectively. The use of Eq. (J3J) is motivated by the semi- 
empirical mass formula for terrestrial nuclei. The bulk 
energy includes not only the volume term but also the 
symmetry term. Since Wb is proportional to the volume, 
we rewrite it as 



W b = w b (u,x in ,n in ,n out )a 3 , 



(4) 



where Wb(u, x ln , n ln , n out ) is the average energy density. 

We assume that the surface energy is proportional to 
the surface area and, hence, depends on the shape of the 
nucleus. Note, on the other hand, that Wb is independent 
of the shape. We rewrite the surface energy in the unit 
cell as 



W s = a{x in ,n" 
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where a(x m , n u 



fc ) is the surface tension, g(u, shape) 



is the relative surface area and shape — SP, C, S, CH, 
SH, G, D, GH or DH. 

The Coulomb energy of the unit cell can be expressed 
similarly as 
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where e is the elementary charge and wc{u, shape) is the 
relative Coulomb energy, which depends on the fraction 
u as well as on the nuclear shape, and is obtained by nu- 
merically solving the Poisson equation for the Coulomb 
potential by a discrete Fourier transform. Since the 
Coulomb energy is proportional to the product of two 
charges divided by their separation and because the to- 
tal charge in the unit cell is proportional to a 3 (Note that 
the charge density is fixed in our model.), the Coulomb 
term, Wc, in the energy expression is proportional to a°. 

Substituting Eqs. g])-© into Eq. we minimize the 
total energy density. This consists of two steps. Firstly, 
minimization with respect to the size of the unit cube, a, 
is performed by a conventional derivation: 
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This result is simply the well-known condition for size 
equilibrium, W s = 2Wc- Eliminating a, we obtain 
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FIG. 2: (a) Relative energy densities, F(u, shape), as a function of u for the periodic bicontinuous morphologies, (b) Differences 
between the relative energy densities of various phases and that of the slab phase, (c) Differences between the average energy 
densities and that of the most stable phase for supernova matter. In all panels, the following notation is used: gyroid (G), 
double-diamond (D), gyroid hole (GH), double-diamond hole (DH), sphere (SP), cylinder (C), slab (S), cylindrical hole (CH) 
and spherical hole (SH). 



Secondly, we minimize Eq. © with respect to the 
shape for a given u. This is performed by simply com- 
paring the energies for different shapes and finding the 
shape that gives the lowest energy. Note that in Eq. ([5]) 
the shape dependence is entirely encapsulated in the rel- 
ative energy density, F(u, shape), denned as 



F(u, shape) = g(u, shape) 2 ^ 3 wc(u, shape 



a/3 



(9) 



As a result, the nuclear shape does not depend on the 
nuclear interaction models, which are encoded in the av- 
erage energy density, 11^(11, x 



n m , n out ), and the surface 



tension, a{x ln ,n 1 ' 



For the conventional pasta phases, the relative energy 
densities are expressed as follows: 
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The numerically determined coefficients, Cbcc = 6.5620 x 
10~ 3 and Chox = 1.2475 x 10~ 3 , are corrections to the 
Wigner-Seitz (WS) approximation, and the subscripts 
"bec" and "hex" represent the body-centered cubic and 
hexagonal lattices, respectively, which are the most sta- 
ble alignments for each nuclear shape 0] • For the hole nu- 
clei, the following relations hold: F(u, CH) = F(l — it, C) 
and F(it,SH) = F(l - u,SP). 

Results. — In Fig. [5[a), we show F(u, shape) as func- 
tions of it for the periodic bicontinuous shapes. It is found 



that the G phase always has a lower energy than the D 
phase. The curves of G (D) and GH (DH) are reflec- 
tions of each other in the vertical line u = 1/2 for the 
following reason. The hole morphology with a fraction 
of u has the same surface area as the normal morphol- 
ogy with a fraction of 1 — it. Hence, relations such as 
g{u, GH) = g(l — u, G) hold. In a similar way, it can be 
easily understood that the hole morphology has a charge 
density and Coulomb potential with opposite signs to 
those of their normal counterparts. Since the Coulomb 
energy is a product of the charge density and Coulomb 
potential, relations such as wc(m,GH) = Wc(l — u, G) 
are satisfied. As a result, the most stable bicontinuous 
phase is G for u < 0.5 and GH for u > 0.5. Note that 
the curves for G and D converge to at the point u = 1, 
which corresponds to uniform matter. In Fig. [2jb), we 
show F(u, shape) — F(u, S) for all the shapes in nuclear 
pasta. We find that the relative energy density of the G 
phase becomes very close to those of the C and S phases 
at the transition point from the C phase to the S phase 
(it = 0.35). The same is true of the transition point from 
the S phase to the CH phase (it = 0.65). 

To evaluate the energy difference quantitatively, we 
must fix the parameters. Bearing in mind the applica- 
tion to the core of supernovae, we adopt a simple es- 
timate using the incompressible liquid drop model em- 
ployed in Ref. [3j. We set n out = and n m = no, where 
no = 0.165 fm~ 3 is the saturation density, and the vol- 
ume fraction is given as it = n/no, where n is the aver- 
age number density of nucleons. The surface tension is 
assumed to be a = 0.73 MeV fm~ 2 , the value that repro- 
duces the properties of isolated finite nuclei in the limit 
of u — * 0. The proton fraction is set to x ln = 0.3. 

The difference between the average energy density for 
each shape and that of the most stable phase, AW/ a 3 , 
for each shape is computed using these parameters and 
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is shown in Fig.[2lc). It can be seen that AW/ a 3 for the 
G phase becomes only ~ 0.2 keV fm -3 at n = 0.35no 
(u = 0.35). This is comparable to the hexagonal lattice 
correction to the WS approximation for cylindrical nuclei 
at u = 0.35. For the matter inside the crust of neutron 
stars, the minimum of AW/a 3 will be even smaller by a 
factor of a few than that for the matter in supernovae. 
Note that the neutron drip is not negligible for matter 
in a neutron star crust and, as a result, the transition 
density will be somewhat larger than n — 0.35no. 

Discussion. — The G phase is not the most stable 
phase for any n (or u) . Considering the tiny difference in 
energy densities, however, it is possible for the G phase 
to exist as a metastable state. In fact, recent studies on 
the dynamics of pasta phases by quantum molecular dy- 
namics (QMD) show that the phase transition from C 
to S is a dynamical process [101 ] and that intermediate 
phases, which are different from any of the known pasta 
phases, may emerge between the C and S phases as well 
as between the S and CH phases These results in- 



dicate the appearance of metastable states, which may 
include the G phase. Since the dynamical stability of 
the G phase cannot be assessed by our model, our argu- 
ment is speculative, but this is an interesting issue to be 
pursued further in the future. 

The coexistence of the G and other phases is even more 
likely at finite temperatures. If this occurs, it should have 
an impact on the thermodynamical properties of nuclear 
matter such as the equation of state. In fact, the mini- 
mum value of AW/a 3 for the G phase is ~ 0.2 keV fm -3 , 
or ~ 3 keV per nucleon, at n — 0.35no, which is much 
smaller than the temperatures of several MeV in super- 
nova cores. Hence, the G phase will almost certainly 
exist as a thermal fluctuation. The surface temperatures 
of some neutron stars are several keV. Since AW/a 3 for 
matter in a neutron star crust is smaller than that for su- 
pernova matter, it is also possible that the G phase will 
appear in neutron star crusts. Note that the shell effect, 
which is also a constituent of the semi-empirical mass for- 
mula but has been estimated to be minor (< 2 keV per 
nucleon) for pasta nuclei [l2j and thus was omitted in our 
model, might also be helpful for the G phase appearance. 

In Ref . [13| , the difference between the energy density 
of the D and S phases was calculated and was found to 
be much larger than that obtained in this paper. This 
discrepancy originates from the fact that in Ref. [3] the 
total energy was not minimized with respect to a. Since 
the energy differences between phases are small for nu- 
clear pasta, the minimization is crucial. 

We have already pointed out an interesting quanti- 
tative similarity between nuclear and polymer systems. 
The G phase in the block copolymer is often observed 
in a narrow range of the volume fraction near ~ 0.35 
This is exactly where F(u, G) becomes very close to 
F(u, C) and F(u, S) in the CLDM. Since Coulomb screen- 
ing cannot be neglected in the polymer system, a sim- 



ple analogy hardly seems applicable. Our results, how- 
ever, appear to suggest some common underlying physi- 
cal principles. [3] 

In this letter we have demonstrated, contrary to con- 
ventional wisdom, that there is a good chance of the G 
phase appearance. There are some issues remaining un- 
resolved in our model, however. Although we assumed 
the shape of the bicontinuous surface a priori, the energy 
density may be further decreased by varying the shape. 
The finite thickness of the surface, which is not taken into 
account in the CLDM, is also expected to be important. 
Therefore, further theoretical investigations with more 
advanced treatments of the surface are necessary. The 
dynamical stability of the G phase, finite-temperature 
corrections to the free energy and the possible implica- 
tions for astrophysical phenomena (e.g., supernova ex- 
plosions, proto-neutron star cooling and pulsar glitches) 
will be also interesting issues. Although there are still 
some controversies, even regarding the appearance of the 
conventional pasta phases, we strongly urge that the new 
type of nuclear pasta is included in these detailed inves- 
tigations. 
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